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Abstract. 

In this study, we examine the sensitivity of long lived tracers to changes in 
the base transport components in our 2-D model. Changes to the strength of the 
residual circulation in the upper troposphere and stratosphere and changes to the low^er 
stratospheric K zz had similar effects in that increasing the transport rates decreased 
the overall stratospheric mean age, and increased the rate of removal of material from 
the stratosphere. Increasing the stratospheric K yy increased the mean age due to 
the greater recycling of air parcels through the middle atmosphere, via the residual 
circulation, before returning to the troposphere. However, increasing K yy along with 
self-consistent increases in the corresponding planetary w^ave drive, which leads to a 
stronger residual circulation, more than compensates for the A^-effect, and produces 
significantly younger ages throughout the stratosphere. Simulations w r ith very small 
tropical stratospheric K yy decreased the globally averaged age of air by as much as 
25% in the middle and upper stratosphere, and resulted in substantially weaker vertical 
age gradients above 20 km in the extratropics. We found only very small stratospheric 
tracer sensitivity to the magnitude of the horizontal mixing across the tropopause, and 
to the strength of the mesospheric gravity wave drag and diffusion used in the model. 
We also investigated the transport influence on chemically active tracers and found a 
strong age-tracer correlation, both in concentration and calculated lifetimes. 

The base model transport gives the most favorable overall comparison w r ith a 
variety of inert tracer observations, and provides a significant improvement over our 
previous 1995 model transport. Moderate changes to the base transport were found 
to provide modest agreement with some of the measurements. Transport scenarios 
with residence times ranging from moderately shorter to slightly longer relative to 
the base case simulated N 2 0 lifetimes that were within the observational estimates 
of Volk et al. [1997], However, only scenarios with rather fast transport rates w'ere 
comparable with the Volk et al. estimates of CFC1 3 lifetimes. This is inconsistent with 
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model-measurement comparisons of mean age in which the base model or slightly slower 
transport rates compared the most favorably with balloon SF 6 data. For all comparisons 
shown, large transport changes away from the base case resulted in simulations that 
were outside the range of measurements, and in many cases, far outside this range. 
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1. Introduction 

Two-dimensional (2-D) chemistry and transport models have been used extensively 
to study the natural and man-induced changes in the middle atmosphere. Of 
fundamental impoitance in 2-D modeling is the proper representation of the zonally 
a\ ci aged transport fields. These include a mean advective component, usually 
represented by the residual circulation of the transformed Eulerian mean (TEM) 
formulation, and a diffusion tensor which parameterizes the eddy flux effects. This 
diffusion is designed to mimic the zonally averaged tracer mixing effects of three- 
dimensional (3-D) planetary wave dissipation, small scale gravity wave breaking, and 
convective overturning. 

Orel tin 1 past two decades, advances in 2-D model transport have involved the 
development of parameterizations of the various eddy processes. For example, the zonal 
mean momentum forcing and horizontal diffusion (7v, /(/ ) induced by planetarv waves can 
be determined empirically, either directly from 3-D meteorological fields, or inferred from 
the diabatic circulation using observed temperatures and heating rates [e.g., Newman 
et al., 1988]. These planetary wave quantities can also be computed theoretically from 
a wave bieaking parameterization [ Garcia . 1991]. Other parameterizations compute 
the wave drag and vertical diffusion (I\ zz ) of breaking gravity w r aves in the mesosphere 
[Lindzen, 1981], and the momentum deposition due to Kelvin wave absorption in 
the tropical upper stratosphere associated with the semiannual oscillation (SAO) 
[Dunkerton, 1979]. Other studies have used long lived tracer observations to estimate 
dynamical quantities in the tropical stratosphere [e.g., Hall and Waugh, 1997; Mote et 
al., 1998], These advances have enabled 2-D models to better simulate the overall zonal 
mean flow structure and general dynamical features seen in observations of long lived 
species such as CH 4 and N 2 0 [e.g., Garcia and Solomon , 1985; Gray and Pyle, 1987; 
Brasseur et al., 1990; Choi and Holton, 1991; Yang et al, 1991; Garcia et al, 1992; Ko 
et al., 1993; Bacmeister et al, 1995; Nightingale et al, 1996; Jackman et al, 1996]. 



5 


Since chemically active species (e.g., CH 4 , N 2 0) can hide some of the transport 
inadequacies in models [e.g., Hall et al., 1999], a more thorough evaluation of model 
transport accuracy can be made by comparing models with observations of inert 
tracers. These include radioactive carbon- 14 ( 14 C) which has a source region in the 
NH middle and high latitude lower stratosphere similar to high speed civil transport 
(HSCT) emissions. Another useful transport diagnostic is the mean age of air (T), which 
corresponds to the mean transit time from the troposphere to a particular location in the 
stratosphere [e.g., Hall and Plumb, 1994], The mean age is related to the stratospheric 
residence time and can be derived from an inert tracer with a quasi-linear increasing 
tropospheric mixing ratio, such as SF 6 [e.g., Elkins et al., 1996; Harnisch et al., 1996]. 
The age of air distribution in 2-D and 3-D stratospheric models has been addressed in 
previous investigations [e.g., Waugh et al., 1997; Bacmeister et al., 1998; Hall et al., 
1999; Li and Waugh, 1999]. 

Recent assessment studies have illustrated the strong dependence of model 
simulations on the stratospheric transport rates used in models [WMO, 1999; Park et 
al., 1999; Kama et al., 1999; IPCC, 1999]. Accurate representation of transport remains 
a major uncertainty in model calculations. In light of this, we recently incorporated 
several of the afore mentioned parameterizations into a comprehensive, empirically- 
based transport formulation in our GSFC 2-D model and obtained significantly better 
agreement with tracer observations compared to previous model versions [Fleming et 
al., 1999]. Subsequently, it is useful to understand how different aspects of the model 
dynamics affect the tracer simulations. 

In the present paper, we build on the work of Bacmeister et al. [1998] and Li and 
Waugh [1999] to further examine the sensitivity of tracer simulations to changes in 
2-D model transport components. We specifically investigate the dependence on the 
residual circulation and vertical diffusion in the upper troposphere and stratosphere, 
the horizontal diffusion and momentum forcing in the stratosphere, and the gravity 
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wave-induced drag and vertical diffusion in the upper stratosphere and mesosphere. 
We will show simulations of various tracers and compare these with observations. This 
allows us to identify the range of model transport rates that are within reasonable 
agreement with the real atmosphere. Simulations of 14 C and age of air are especially 
useful in diagnosing model accuracy in transporting and accumulating HSCT exhaust in 
the stratosphere. This will provide a basis for the companion paper which explores how 
the variations in the model transport affects the simulations of an HSCT perturbation. 

2. GSFC-2D Model 

The 2-D model at NASA/Goddard Space Flight Center (GSFC) was originally 
described in Douglass et al. [1989] and Jackman et al. [1990]. Recent improvements to 
the model have been discussed in Jackman et al. [1996], and upgrades to the transport 
formulation are described in Fleming et al. [1999]. Details of this new transport 
formulation are contained in the Appendix. We note that the bottom boundary 
condition used to compute the model residual circulation (equation A3 in the Appendix) 
is now more physically realistic than that used in recent assessment studies. This 
resulted in a somewhat slower stratospheric circulation and an older mean age of air 
than those published in the recent assessment reports [Kama et al., , 1999; Park et al., 
1999; Hall et al., 1999]. 

We have updated the model reaction rates and photolysis cross sections to the 
latest Jet Propulsion Laboratory (JPL) 1997 recommendations [DeMore et al., 1997]. 
We have also made some minor modifications to our model chemistry scheme. H 2 0 2 
is now transported, and the HO x species (H -I- OH + H0 2 ) are now computed and 
transported using the family approach following Douglass et al. [1989]. The family 
members are computed from photochemical equilibrium and normalized to the total 
family concentration. In a similar fashion, the hydrocarbon products of CH4 oxidation 
are now also computed as a family (CH X = CH 3 + CH 3 0 2 + CH 3 0 + CH 2 0 + CH 3 OOH). 
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This replaces our previous methodology in which the individual HO x and CH X species 
were computed separately assuming photochemical equilibrium. This new formulation 
ensures that the total hydrogen and carbon atoms are conserved throughout the model 
run. However, we note that these changes did not impact on the model simulations 
shown in this study. 

3. Results 

We now show model simulations of several dynamical tracers from a variety of 
transport sensitivity tests, listed in Table 1. As discussed in Fleming et al. [1999], 
simulations with the base line model transport give the best overall comparisons with 
tracer observations. In the experiments shown here, we will illustrate how moderate and 
large changes away from this baseline transport affect the model simulations of long 
lived tracers. We will discuss how changes in the circulation, K yy , and K zz fields affect 
the simulations separately, and then in various combinations. We will also compare the 
results from our latest 1999 model with those obtained using the previous 1995 version 
of our model transport. 

3.1. Sensitivity to the circulation 

We first examine the model sensitivity to changing the residual circulation 
throughout the troposphere and lower stratosphere (the K yy and K zz fields remain 
fixed). In these experiments, we have increased and then decreased the circulation 
strength away from the baseline by moderate and then large amounts in the region 
below 35 km. These are defined as scenarios B, C, D, and E as listed in Table 1. 

Figure 1 shows observations of 14 C at 31°N for July 1964 and July 1966, along 
with time dependent model simulations for the base transport (scenario A) and the 
four circulation sensitivity scenarios. The profiles for July 1964 correspond to 9 months 
after initialization. Both the observed and simulated peaks in 14 C at 20-25 km are 
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still fairly robust, with a strong decrease with height into the middle stratosphere. 
The weakest circulation gives the largest 14 C concentrations below the peak, with the 
strongest circulation producing the largest 14 C above 28 km during this time. Here the 
concentrations are directly related to the strength of the circulation which transports 
14 C up from the source region near 20 km. Also, all of the circulation scenarios simulate 
larger than observed 14 C concentrations above ~26 km during this time. By July 1966 
(2 years and 9 months after inialization), the maximum has been dissipated with very 
little vertical gradient seen above 22 km in the simulations and observations. At this 
time, removal of 14 C from the stratosphere is proportional to the circulation strength so 
that more 14 C is maintained in the stratosphere with a weaker circulation, similar to the 
findings of Jackman et al. [1991]. Qualitatively similar results are seen in simulations of 
strontium 90 (not shown) in the tropics and SH midlatitudes. 

The annually averaged mean age of air derived from SF6 for the baseline and 4 
circulation scenarios (B, C, D, and E) is shown in Figure 2. In this study we use the 
global mean SF6 time series at the surface [ Geller et al., 1997] as the reference for 
computing T in the stratosphere. The shapes of the contours in Figure 2 are typical 
of a long lived tracer illustrating the Brewer-Dobson circulation, with rising motion 
in the tropical stratosphere and descent in the extratropics. The age throughout 
the stratosphere becomes progressively older with a weaker circulation, although the 
general shapes of the age contours do not change appreciably with circulation strength. 
Ages of greater than 9 years are simulated in the upper stratosphere with the weakest 
circulation. This age dependence on the circulation strength is also illustrated in Table 1 
which lists the annual/global average T at 19, 27, and 35 km for each transport scenario. 

Figure 3 shows the mean age of air at 20 km for Oct/Nov 1994 derived from 
ER-2 SF 6 data taken from the Airborne Chromatograph for Atmospheric Trace Species 
(AC ATS) instrument [Elkins et al., 1996], along with the different model circulation 
scenarios. Figure 4 shows vertical profiles of the mean age at low, middle, and high 
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latitudes from the model simulations and ballon SF 6 observations of Harnisch et al. 
[1996] (asterisks) and the Observations of the Middle Stratosphere (OMS) campaign 
Lightweight Airborne Chromatography Experiment (LACE) instrument (F.L. Moore 
et al.. manuscript in preparation. 1999) (triangles). Note that mean ages derived from 
C'Oj measurements [e.g.. Andrew . s el al.. 1999] were found to be similar to those derived 
from SF (i data al 20 km [Hall it al.. 1999]. 

Figures 3 and -1 again show that increasing the circulation strength decreases the 
stratospheric T everywhere above* 15 km. The age differences between the different 
circulations are largest at the highest altitudes, with a ~-6 year difference at 33 km 
between the very weak and very strong circulations. The overall vertical gradients 
and the latitudinal gradients between the tropics and mid-latitude.' also decrease with 
increasing circulation strength 1 his shows lliat such gradients ean be modified by 
the residual circulation its well as the horizontal diffusion rales. Note also that the 
dependence of the vertical age gradients on the circulation and the overall vertical 
structure of the midlalitude age distribution is very similar to that seen in the July 1966 
U C simulation at 31°X in Figure 1. 

Figures 1-4 illustrate the range of results that are obtained with the various 
circulation strengths, and how these compare against observations. While the baseline 
transport provides the best overall model-data comparison, both the moderately weak 
and moderately strong circulations (scenarios C and D) are within reasonable agreement 
with the observations in some of the comparisons. For example at 20 km in the 
tropics, the base and strong circulations compare well with the data, whereas the weak 
circulation indicates ages that are too old. However, the weak circulation is within 
the range of observations at 65°N in Figure 4, and provides the best agreement above 
23 km with the OMS midlatitude data (triangles, middle panel). However, the base 
model agrees best with the midlatitude data of Harnisch et al. [1996] (asterisks). We 
note that the two midlatitude data sources show age differences of about 1 year in the 
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middle stratosphere, even though they are for similar latitudes and seasons. Some of 
this difference may be due interannual dynamical variability since the data were taken 
during different years (1996 versus 1993). Such variability is not accounted for in the 
model simulations. 

It appears that the very weak and very strong circulations (scenarios B and 
E) produce ages that are outside the range of observations, and in some cases, far 
outside this range. An exception to this is the 65°N profile in which there are a few 
measurements of very old air (7-9 years) at 20-25 km thought to be remnants of the 
polar vortex [ Ray et al., 1999]. Such phenomena cannot be resolved by a zonal mean 
model, except with an artificially very weak circulation as shown in Figure 4. We 
note that it is possible that such observations of old ages are an artifact, in part, of a 
mesospheric photochemical loss of SF 6 [e.g., Hall and Waugh , 1998] which would be 
reflected in polar vortex air. 

As a final transport diagnostic, Figure 5 shows vertical profiles of C0 2 plotted 
against N 2 0 in the SH polar vortex region from the ER-2 measurements [Boering et 
al., 1996; Andrews et al, 1999] and the model simulations. The C0 2 -N 2 0 relationship 
in the stratosphere has been discussed previously [Boering et al., 1996; Strahan et al. , 
1998]. The tropospheric source of C0 2 has a secular trend along with a strong seasonal 
cycle in the NH. The source of N 2 0 is seasonally invariant and undergoes only a small 
secular increase. At high southern latitudes where the C0 2 seasonal amplitude in the 
lower stratosphere is quite small, it would be expected that variations in C0 2 for a given 
value of N 2 0 will be due, in part, to the transport of the differing trends of the C0 2 and 
N 2 0 sources. Also, N 2 0 undergoes photochemical loss in the upper stratosphere which 
partially offsets the N 2 0 increase caused by a stronger circulation and larger mass flux 
from the tropical troposphere. C0 2 has no photochemical sink in the stratosphere so 
that it will be much more sensitive to transport changes compared to N 2 0. 

Figure 5 is representative of most latitudes and seasons. This shows that the 
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CO 2 — N 2 O correlation is strongly sensitive to the circulation such that for a given value 
of N 2 0, CO 2 increases with an increasingly stronger circulation. With a very weak 
circulation, the vertical (and horizontal) gradients in CO 2 are significantly stronger 
relative to the baseline case. These gradients become w r eaker with increasing circulation 
so that with the very strong circulation, CO 2 is transported rapidly enough throughout 
the stratosphere to produce little or no horizontal or vertical gradients. Figure 5 
illustrates that the base transport shows the best overall consistency with the data, with 
the other circulations being fairly far outside of the range of data, at least above N 2 0 
levels of 275 ppbv. Below this level, the strong and very strong circulations (scenarios 
D and E) provide a better comparison with the data. However, incomplete resolution 
of processes near the tropopause in the model and/or localized effects in the data could 
also affect the model-measurement comparison in this region. 

Given that CO 2 acts as a conserved tracer in the stratosphere and has a linear 
increase at the ground, Figure 5 can also be thought of as a vertical profile of the mean 
age. Since the seasonal cycle amplitude in C0 2 is quite small at high SH latitudes, the 
profiles in Figure 5 reflect the fact that the linear trend in the CO 2 source is transported 
through the stratosphere more rapidly with a stronger circulation. For example at the 
lowest N 2 O values of 125 ppbv (corresponding to 22-25 km), the very weak and very 
strong circulations show a CO 2 difference of 8.6 ppmv. Assuming a linear increase in 
CO 2 of 1.4-1. 5 ppmv/yr at the ground, this corresponds to an age difference of 5. 7-6.1 
years between the very weak and very strong circulations. This is consistent with the 
differences in T derived from SFg in the SH polar lower stratosphere shown in Figure 2. 

3.2. Sensitivity to lower stratospheric K zz 

Observational studies [e.g., Hall and Waugh, 1997; Mote et al., 1998], have shown 
that the vertical diffusion ( K zz ) in the lower tropical stratosphere is .01-. 02 m 2 /s, and 
we have adopted this value in our base model simulation. In this section, we investigate 
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the sensitivity to changing the K zz values in the lower stratosphere away from this 
baseline value. 

Figure 6 shows the model simulations of 14 C in July 1966 for minimum lower 
stratospheric K zz values of 0.001, 0.1, and 1 m 2 /s (scenarios F, G, and H). Also plotted 
here are the base model scenario A and observations identical to Figure 1. We found 
that using a minimum stratospheric K zz smaller than 0.001 m 2 /s does not change the 
model response in the diagnostic tracers so that this value appears to be the lower limit 
of K zz in our model. These scenarios are consistent with the circulation changes in 
Figure 1 in that a slower (faster) transport maintains more (less) 14 C throughout most 
of the stratosphere. Values of K zz between 0.001 and 0.1 m 2 /s cause relatively small 
changes in the 14 C simulation above 18 km, with the baseline or K zz value of 0.1 m 2 /s 
providing reasonable agreement with the data. There is very little difference among 
scenarios A, F, and G below 18 km, with all 3 providing good agreement with the data. 
Using a K zz of 1 m 2 /s (scenario H) causes mixing within and removal of material from 
the stratosphere that is much too rapid. This results in simulations that are far outside 
the range of observations as seen in Figure 6. 

The annually averaged mean age distributions for the baseline and K zz scenarios are 
shown Figure 7. Figures 8 and 9 show the latitudinal and vertical profiles of mean age 
for these simulations, along with the base case and observations as in Figures 3 and 4. 
Similar to the circulation changes, T throughout the stratosphere becomes progressively 
younger with increasing K zz . Using a minimum K zz beween 0.001 and 0.1 m 2 /s varies 
the age of air by at most 1 year, but the resulting simulations are within the range 
of observations. An exception to this is the vertical profile at midlatitudes where the 
smallest K zz agrees quite well with the Harnisch et al. [1996] profile, but is still ~ 1 year 
too young above 23 km compared to the OMS data. Similar to the circulation changes, 
the latitudinal and vertical age gradients decrease with increasing K zz . 

Consistent with the 14 C simulations, using the very large K zz of 1 m 2 /s (scenario H) 
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simulates ages that are much too young everywhere, and gives latitudinal and vertical 
gradients in V that are much weaker than observed. The simulations with this large 
K zz are consistent with the model-measurement intercomparison analysis of Hall et al. 
[1999] who concluded that a K zz of 1 m 2 /s is unrealistic and that 0.1 m 2 /s is close to 
the expected upper limit on K zz . We note also that this tracer dependence on K zz is 
similarly seen in the simulated C0 2 - N 2 0 relationship (not shown). 

3.3. Sensitivity to mesospheric gravity wave effects 

The model gravity wave parameterization computes diffusion, used as the model 
mesospheric K zz field, and wave drag which has a direct effect on the model residual 
circulation (see Appendix). The relative importance of vertical diffusion and vertical 
advection by the residual circulation in controlling the water vapor distribution in the 
mesosphere has been debated previously [e.g., Garcia and Solomon, 1985; Holton and 
Schoeberl, 1988; Smith and Brasseur, 1991; Nedoluha et al., 1996]. Figure 10 shows the 
seasonal variation in H 2 0 at 80 km as measured by the Halogen Occultation Experiment 
(HALOE - version 19) instrument onboard the the Upper Atmosphere Research Satellite 
(UARS). In our model, vertical advection is the dominant process controlling H 2 0 
throughout most of the mesosphere, and is largely responsible for the H 2 0 seasonal 
cycle [Chandra et al., 1997]. Figure 10 shows that the base model scenario resolves the 
overall seasonal cycle seen in the HALOE data reasonably well. 

Vertical diffusion, which maximizes during the solstice seasons when the westerly 
and easterly midlatitude jets are strongest, plays a secondary role in determining the 
extratropical H 2 0 seasonal variations in our model. This diffusion acts to increase 
(decrease) H 2 0 in the upper (lower) mesosphere. In the upper mesosphere, the diffusive 
transport therefore enhances the advectively-driven water vapor increase in summer. 
Conversely, diffusion partially offsets the advectively-induced H 2 0 decrease in winter. 
These effects maximize at 80 km, resulting in secondary maxima during mid-winter at 
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midlatitudes. This effect is somewhat more pronounced in the SH winter because the 
stronger zonal westerlies, and hence the larger difference between the background zonal 
w T ind and the gravity wave phase speed, allow for larger diffusion compared to the NH 
winter. These secondary mid-winter maxima are seen in the HALOE data in Figure 10, 
and it appears that our model captures these features reasonably w r ell, including the 
larger secondary maximum in the SH. We note that horizontal advection by the residual 
circulation has a substantially smaller influence on mesospheric H 2 0 compared to either 
the vertical advection or vertical diffusion. 

To test how these processes affect tracers in the stratosphere, we ran two additional 
scenarios with weak and strong mesospheric gravity wave-induced drag, which affects 
the circulation, and diffusion, relative to the base case (scenarios I and J). As seen in 
Figure 10, these scenarios have a large effect on mesospheric H 2 0. Stronger gravity wave 
effects lead to larger water concentrations at 80 km due to the net effect of increasing 
the upward transport of H 2 0 from the stratosphere. The base model gives the most 
favorable overall comparison with the HALOE data, especially at midlatitudes. The 
strong gravity wave scenario (J) also shows some reasonable agreement, especially at 
the equator. However, the weak gravity wave case (I) simulates very small mesospheric 
H 2 0, especially in the winter extratropics, and compares rather poorly with the HALOE 
data. 

Despite the large induced changes in the mesosphere, these scenarios exert very 
little influence on tracer distributions in the lower and middle stratosphere in our model. 
The changes in global and annual mean age below 35 km between scenarios I and J and 
the base case were at most 4-5 months (Table 1). Furthermore, only very small changes 
were seen in the model simulations of 14 C below 35 km. 

3.4. Sensitivity to stratospheric K yy 

Figure 11 shows distributions of mean age for several different K yy scenarios 
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along with the baseline model. Here we have decreased and increased the K yy values 
everywhere above the tropopause by constant factors of 2 and 5 (scenarios K, L, M, 
and N in Table 1). As expected, the latitudinal gradients in T become progressively 
weaker with increasing horizontal diffusion. It is also apparent that the mean age 
throughout the stratosphere above 16 km becomes increasingly older with increasing 
Kyy (see also Table 1). This is due to the fact that with larger K yy , air parcels undergo 
more recycling through the middle atmosphere via ascent in the tropics and descent in 
the extratropics, before returning to the troposphere, and hence achieve an older mean 
age [e.g., Bacmeister et al., 1998; Hall et al., 1999], Ages greater than 12 years are 
simulated in the upper stratosphere for the largest K yy case (scenario N, Figure lib). 

The rate of horizontal mixing between the tropics and midlatitudes has been 
examined in previous studies [e.g., Plumb, 1996]. Our base model tropical-subtropical 
K yy field varies spatially and seasonally, ranging from 10 8 to 10 10 cm 2 /s (see Appendix). 
The very large K yy scenario (N) includes large K yy throughout the tropics and is akin 
to the “global diffuser” model discussed by Plumb, [1996]. To further test the tracer 
sensitivity to our model horizontal diffusion across the tropics, we ran an additional 
scenario (Q) in which K yy was set to a very small value of 10 8 cm 2 /s throughout the 
tropical stratosphere (15°S-15°N, 15-50 km) for all seasons. This is similar to the 
“tropical pipe” model of Plumb, [1996]. The resulting mean age contours in Figure Ilf 
show very large horizontal gradients between the tropics and midlatitudes reflecting a 
very isolated tropical region. T is younger throughout the stratosphere in this case, with 
differences of almost 2 years in the high latitude middle stratosphere relative to the base 
case. Note also the small isolated region of T > 5 years in the very lower stratosphere 
over the Antarctic. 

Latitudinal variations of T for these K yy scenarios are shown in Figure 12 which 
again illustrates the general increase in mean age with increasing K yy . At this level, 
it appears that the simulation is in reasonable agreement with the observations for 
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moderate changes in eddy diffusion. Even the smallest K yy case (scenario K) shows 
some agreement with the data at midlatitudes. However, the largest K yy case simulated 
T that is much older than observed at almost all latitudes. At 20 km, the small tropical 
K yy scenario (Q) simulates ages in the tropics that are almost 1 year younger relative to 
the base case, due to the lack of in-mixing of older midlatitude air. The corresponding 
subtropical gradient is much stronger in this case, with slightly younger ages again 
simulated at high latitudes relative to the base case. 

The vertical profiles in Figure 13 show the strong K yy - age dependence above 20 km 
which is similarly seen in C0 2 - N 2 0 profiles (not shown). The extreme scenarios (K 
and N) appear to be far outside the range of data above about 25 km. At midlatitudes, 
the base model agrees with the Harnisch et al., data, while the moderately large case 
(scenario M) agrees best with OMS data. A similar finding was seen in the circulation 
tests in section 3.1. The vertical gradient in T above 20 km also shows a strong 
dependence on horizontal diffusion such with the gradient weakens with a progressively 
smaller global K yy . At middle and high (not shown) latitudes, the vertical gradient 
is very close to zero for the smallest K yy case. This illustrates the dominance of the 
residual circulation, which is primarily downward in this region, in controlling the long 
lived tracer distributions in this case. We note that F from the small tropical K yy 
scenario (Q) (not shown in Figure 13) is nearly identical to the K yy /2 case (scenario L) 
at these latitudes. As the recycling of air throughout the stratosphere is inhibited in the 
reduced tropical K yy case, F at all latitudes above 20 km is decreased. Mean ages are 
younger by as much as 1.5 years at 35 km at 65°N relative to the base case. 

Thus far, we have treated changes in K yy separately from other changes in the 
model transport. However in reality, changes in K yy are due to changes in wave driving 
which simultaneously affect the residual circulation. To investigate this coupling, we ran 
an additional scenario in which K yy has been decreased by a factor of 5 as in scenario 
K, along with a corresponding factor of 5 decrease in the wave drive (also above the 
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tropopause), which weakens the circulation. We also ran a case in which both K yy and 
wave drive were increased by a factor of 5. These are listed as scenarios O and P in 
Table 1. 

The age of air response to increasing or decreasing K yy and wave drive self- 
consistently is illustrated in Figures llg-h. In these cases, the shapes of the contours 
are similar to those in which only K yy has been changed (Figures lla-b). However, 
additionally reducing the wave driven circulation increases the age everywhere, and 
produces isolated regions of old air in the high latitude lower stratosphere in both 
hemispheres. Conversely, increasing the circulation together with the K yy field in Figure 
llh drastically reduces the mean age throughout the stratosphere. 

Latitudinal profiles of these scenarios at 20 km are shown in Figure 14. Again, 
decreasing only K yy (dashed-triple dotted line identical to Figure 12) results in younger 
ages everywhere compared to the base case. As discussed in section 3.1, slowing the 
circulation produces older ages throughout the stratosphere. The net effect of including 
both the weaker wave driven circulation and K yy (dashed-single dotted line in Figure 
14) produces younger air in the tropics, older air in the extratropics, and a stronger 
subtropical age gradient, relative to the base case. 

Figure 14 also shows the dramatic increase in T throughout the stratosphere 
obtained by increasing K yy (dotted line, identical to Figure 12). However, increasing 
both Kyy and the wave drive (dashed line) has the opposite effect, with the mean ages 
being much younger everywhere and the age gradients being much weaker compared to 
either the base scenario or the observations. Note that because of the small inherent 
wave drive in the tropical stratosphere, decreasing this quantity in the tropics to be 
dynamically consistent with the very small K yy used in scenario Q gave virtually the 
same mean age response compared to decreasing K yy only. 

We note that the combined K yy and wave drive changes shown here are a somewhat 
incomplete representation of the true atmosphere. In reality, changes in K yy and the 
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wave driven circulation would also change T, the diabatic heating rates, and the u 
structure, which in turn will affect the wave propagation charcteristics responsible for 
the changes in K yy and wave drive. Quantifying such changes in a completely consistent 
manner is difficult in our model given the empirical nature of the transport fields. To 
be fully self-consistent would require using an interactive chemical-radiative-dynamical 
2-D model, such as reported in Bacmeister et al. [1998]. 

As a final K yy sensitivity test, we investigated the model response to changes in 
horizontal diffusion across the tropopause. Recent reports have noted the importance 
of cross-tropopause transport in model simulations of HSCT emissions [e.g., Park et 
al., 1999]. It was suggested that 2-D and 3-D models simulate a different correlation 
between mean age and NO y accumulation, perhaps because of the inherent differences 
between 2-D and 3-D model formulations of cross-tropopause transport, which can 
have large zonal asymmetries. In our 2-D model, the tropopause, which is a function 
of latitude and season and is based on NCEP data, separates very fast and very slow 
vertical mixing ( K zz ) regimes. There is also cross-tropopause transport horizontally 
via the K yy field where the tropopause height changes with latitude. This occurs 
primarily at midlatitudes between ±25° — 60°. To test how this horizontal transport 
affects the diagnostic tracer simulations, we ran an additional model scenario (R) in 
which the cross-tropopause K yy values were significantly increased over the baseline 
case (the K zz values remain unchanged). In this case, simulations of long lived tracers 
were nearly identical to the baseline case. Mean ages were at most 2 months older 
throughout the stratosphere in scenario R compared to the base case (Table 1). As 
shown in the companion paper, the T - NO y correlation also changed only very slightly 
in this case. This illustrates the very small tracer sensitivity to the rate of horizontal 
cross-tropopause transport in our model simulations. 
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3.5. Sensitivity to combined transport changes 

In this section we explore how simultaneous changes to the circulation, K yy , and 
K zz fields affect the long lived tracer simulations. Scenario S combines transport changes 
which all contribute to a moderately long stratospheric residence time, as reflected in 
the mean age in Figures 15a, 16 and 17 (see also Table 1). These include the weak 
circulation (scenario C), small lower stratospheric K zz (scenario F), and moderately 
large K yy (scenario M). This scenario gives very old ages of greater than 12 years 
in the upper stratosphere and mesosphere, with a strong vertical age gradient in the 
lower stratosphere (Figure 17, dashed-triple dotted line). Although the latitudinal age 
gradients in the middle and upper stratosphere are generally weaker than the base case 
due to the larger K yy , the gradients are actually a bit stronger in scenario S in the lower 
stratosphere (Figure 16). 

Scenario T combines the transport changes which contribute to a moderately 
short stratospheric residence time. These include the strong circulation (scenario D), 
moderately large K zz (scenario G), and moderately small K yy (scenario L). The young 
ages produced by this scenario are seen in Figure 15b, along with weaker latitudinal 
gradients in the lower stratosphere (Figure 16, dotted line) and weaker vertical gradients 
throughout the stratosphere (Figure 17, dotted line). Figures 16 and 17 show that at 
levels above the very lower stratosphere, both scenarios S and T give tracer simulations 
that are significantly different from the observations. 

To test the extreme limits of these sensitivity tests, we combined the transport 
scenarios w’hich had the very longest and very shortest stratospheric residence times 
discussed in this study. Scenario U combines the very weak circulation (scenario B), 
small K zz (scenario F), and very large K yy (scenario N), with the opposite combination 
(scenarios E, H, and K) used for scenario V. As shown in Figures 15c-d and in Table 1, 
mean ages in the middle stratosphere are > 20 years for scenario U and < 2.5 years for 
scenario V . These simulations are even further outside the range of measurements than 
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those given by scenarios S and T in Figures 16-17. These cases represent the absolute 
extremes of stratospheric mean age and residence time that we obtained in this study 
by varying the advective and diffusive transport rates in our model. 

3.6. 1995 model transport 

As a final scenario, we ran the tracer imulations using the previous 1995 version of 
our model transport (scenario W). This version of the model was used in the 1995 NASA 
HSCT assessment [Stolarski et al., 1995] As seen in Figures 16-17, the new 1999 model 
transport (base scenario A) simulates tracers that are in significantly better agreement 
with observations compared to the previous 1995 transport [see also Fleming et al., 
1999]. The contour shapes of F are also significantly different between the 1999 and 1995 
models (Figures 15e-f), with the 1995 model simulating greater asymmetry between 
the hemispheres and weaker gradients overall than the 1999 model. Furthermore, the 
gradients and absolute values of 14 C and mean age in the 1995 model are significantly 
underestimated, relative to observations, indicating that the transport rates are too 
fast in this older model version. The simulated annual/global mean ages are 40-60% 
younger throughout the stratosphere compared with the 1999 model transport. The 
1995 model simulations in Figures 16-17 are actually similar to those of the combination 
short residence time scenario (T). 

There are substantial differences between our 1999 and 1995 model transport 
formulations which are reflected in these simulations. While there have been changes 
to the mean circulation, much of the differences are due to changes in the meridional 
and vertical diffusive transport. The 1999 model is generally less diffusive in a globally 
averaged sense, and exhibits greater isolation between the troposphere and stratosphere, 
i.e., a sharper change in K zz across the tropopause. There is also a sharper transition in 
Kyy, and hence greater isolation, between the tropics, subtropics, and midlatitudes in 
the stratosphere in the 1999 model. 
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3.7. Chemically active tracers 

Previous studies have examined how transport variations affect tracers that are 
chemically active [e.g., Hall et al, 1999; Schoeberl et al, 1999; Li and Waugh, 1999]. 
Because of the rapid photochemical time scales in the upper stratosphere, constituents 
such as N 2 O and Cl y are less dependent on transport rates compared to inert tracers 
such as SF 6 . Therefore, chemically active tracers can hide model transport changes and 
hence, transport deficiencies. Figure 18 shows the global and annually averaged Cl y and 
N 2 O concentrations at 27 km plotted against mean age for the 23 transport scenarios 
discussed in this study. This level is characteristic of other stratospheric altitudes 
and shows a strong correlation of mean age with both Cl y and N 2 0. Faster transport 
rates result in younger ages, and yield smaller concentrations of Cl y (less conversion of 
inorganic chlorine to Cl y ), and larger N 2 O values (less photochemical destruction). The 
base scenario (depicted by the large letter “A” in Figure 18) appears near the middle of 
the different transport scenarios. 

Figure 18 also shows that the correlation pattern is not linear, but rather shows 
an asymptotic behavior. This is consistent with previous studies [e.g., Hall et al, 1999; 
Schoeberl et al, 1999; Li and Waugh, 1999], suggesting that very slow transport rates 
which give long stratospheric residence times strongly affects inert tracers such as F. 
However, photochemically active species such as N 2 O and Cl y show progressively less 
sensitivity to decreasing transport rates, since most if not all of the N 2 0 has been 
photochemically destroyed, and the inorganic chlorine converted to Cl y in very old air. 

The response of chemically active tracers to the stratospheric residence time, or 
lifetime, is also seen in Figure 19. Here we plot the correlation of mean age with 
the stratospheric steady state lifetimes of CFC1 3 and N 2 0 for the different transport 
scenarios (see also Table 1). Recently, Volk et al. [1997] estimated the lifetimes of several 
source gases using ER-2 aircraft observations, and found significantly shorter lifetimes 
compared to those derived from 2-D model calculations [e.g., WMO , 1995]. The Volk et 



22 


al. [1997] best estimate of the lifetimes of CFCI3 and N2O, along with the 1 standard 
deviation uncertainties, are plotted with our model results in Figure 19. 

As in Figure 18, the model simulations show a strong correlation between mean 
age and the CFC1 3 and N 2 0 lifetimes. With faster transport rates there is a larger 
constituent concentration in the stratosphere, but a greater increase in the constituent 
loss rates compared to the base case, so that the net effect is to decrease the lifetime. 
The correlation with mean age is almost completely linear for these tracers. Our base 
scenario reveals a stratospherically averaged (pressure weighted) mean age of 2.9 years, 
with lifetimes for CFCI3 and N 2 0 of 68 and 136 years, respectively. Scenario V simulates 
the shortest stratospheric residence time and youngest mean age, with lifetimes of 77 
years for N 2 0 and 27 years for CFCI3. Conversely, the longest residence time and oldest 
T is simulated by scenario U, which gives lifetimes for N 2 0 and CFCI3 of 461 years and 
231 years, respectively. 

The majority of our model transport scenarios simulate N 2 0 lifetimes that are 
within the range of the Volk et al. [1997] estimates (98-146 years). These include several 
scenarios with very fast transport rates, e.g.^ the very strong circulation (E) and fast 
stratospheric vertical mixing (H), even though the mean age and 14 C distributions for 
these scenarios were far outside the observational range. The base scenario N 2 0 lifetime 
is somewhat longer than the best empirical estimate, but is still within the range of 
observations. For the CFCI3 lifetimes, only the rather short residence time scenarios are 
within the range of the Volk et al. estimates, and the base scenario CFCI3 lifetime of 68 
years is ~50% longer than the observations. Scenarios with very long residence times, 
e.g. the very weak circulation (B), fast horizontal mixing (N), and the combination 
scenarios S and U, revealed lifetimes for both N 2 0 and CFCI3 that are much longer 
than observed. Interestingly, the 1995 model transport (scenario W) appears to give 
the best overall agreement with the observations for both the N 2 0 and CFCI3 lifetimes. 
However as seen in Figures 16-17, this scenario gave significantly poorer agreement with 
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other tracer observations compared to the base 1999 transport scenario (A). 

4. Conclusions 

In this study, we have examined the effects on stratospheric long lived tracers of 
relatively moderate and large changes to the base transport fields in our 2-D model. 
Changes to the strength of the residual circulation in the upper troposphere and 
stratosphere had an effect similar to changing the lower stratospheric K zz values: 
increasing the transport rates decreased the overall stratospheric residence time, and 
hence the mean age, and increased the rate of removal of material from the stratosphere. 

Increasing the stratospheric horizontal diffusion decreased the latitudinal tracer 
gradients, but increased the overall stratospheric mean age. This is due to the fact that 
with a larger K yy , air parcels undergo more recycling through the middle atmosphere 
via ascent in the tropics and descent in the extratropics, before returning to the 
troposphere. Simulations with very small tropical stratospheric K yy resulted in a very 
isolated tropical region and reduced the mean age throughout the stratosphere, with T 
younger by as much as 25% in the mid-upper stratosphere. Also, large scale reductions 
in Kyy strongly reduced the vertical gradients in T above 20 km in the extratropics, 
illustrating the dominance of the residual circulation in controlling the long lived tracer 
distributions with relatively small global K yy values. 

We found that significant variations in the strength of the mesospheric gravity wave 
drag and diffusion strongly influenced the simulated seasonal cycle of water vapor in the 
upper mesosphere, but caused only very minor changes to 14 C and the mean age below 
35 km. Similarly, the model simulations in the stratosphere showed only very small 
sensitivity to the magnitude of the horizontal diffusive transport across the tropopause. 

Incorporating K yy changes along with self-consistent changes in the corresponding 
planetary wave drive, and hence, the residual circulation, had a cancellation effect, and 
produced significantly different tracer simulations compared to changing only the K yy 
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rates. While increasing K yy globally increases the mean age, a stronger wave drive, 
and hence, an enhanced circulation implied by a larger A yy acts to decrease the mean 
age, with the net effect giving significantly younger ages everywhere relative to the base 
case. The shortest stratospheric residence time and youngest mean age obtained in this 
study were simulated by combining a strong circulation, large K zz , and small K yy field. 
Conversely, combining a weak circulation, small lower stratospheric K zz , and a large 
stratospheric K yy field yielded the longest stratospheric residence time and oldest mean 
age. The upper and lower limits of the annual/globally averaged middle stratospheric F 
obtained in this study ranged from ~1 year to slightly more than 20 years. 

The model simulations revealed a strong correlation between mean age and 
chemically active tracers, both in concentration and calculated lifetimes. The correlation 
between T and both N 2 O and Cl v concentrations showed an asymptotic behavior. This 
is indicative of the fact that as transport rates decrease and stratospheric residence 
times and mean age increase, N 2 O arid Cl y , which are photochemicallv active, become 
progressively less sensitive to transport. 

The base model transport provided the most favorable overall comparison with a 
variety of observations, and represents a significant improvement over our 1995 transport 
formulation. Moderate changes to the base transport were found to provide modest 
agreement with some of the measurements. Scenarios with residence times ranging 
from moderately shorter to slightly longer relative to the base case simulated N 2 O 
lifetimes that were within the observational estimates of Volk et al [1997]. However, 
only scenarios with rather fast transport rates were comparable with the Volk et al 
estimates of CFCI3 lifetimes. This is inconsistent with model-measurement comparisons 
of mean age in which the base model or scenarios with transport rates slightly slower 
than the base case compared the most favorably with balloon SF 6 data. We also found 
that large transport changes away from the base case resulted in simulations that were 
outside the range of measurements, and in many cases, far outside this range. 
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These transport scenarios can be used to explore the sensitivity to dynamics of 
model simulations of an HSCT perturbation. This will be addressed in part 2 of this 
study. 
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Appendix: GSFC 2-D Model Transport Description 

Earlier versions of the model transport derived the residual circulation (v*,w*) 
using climatological temperatures and heating rates in the TEM thermodynamic 
and continuity equations [e.g., Dunkerton, 1978]. The planetary wave-induced K yy 
distribution was determined self-consistently from the circulation using the residual 
needed to balance the zonal momentum equation [Newman et al., 1988]. A simple, 
weakly height-varying K zz field with no latitudinal or seasonal variations was used to 
describe turbulent diffusion induced by gravity wave breaking in the mesosphere. These 
dynamical fields were able to capture some general features seen in long-lived tracer 
observations. However, this transport did not explicitly include planetary or gravity 
wave effects. Therefore, many details of the observed zonal mean transport were not 
resolved in the model simulated tracer fields. 

We have upgraded the model transport formulation which now generally follows 
that of interactive 2-D models [Garcia and Solomon , 1983]. However, the new dynamical 
fields are based on empirical data sets as opposed to being computed interactively in 
the model. This scheme uses explicit diagnostic calculations of planetary and synoptic 
scale wave effects, along with parameterizations to account for gravity wave and Kelvin 
wave effects in which the zonal mean flow is constrained to observations. Compared to 
our previous methodologies, this new formulation provides an improved representation 
of the zonal mean dynamics [Fleming et al., 1999]. As with previous model versions, 
the transport is climatological and does not account for the quasi-biennial oscillation or 
other interannual variability. Details of the new formulation are described below. 

A.l. Zonal Mean Formulation 

Following the methodology originally outlined in Garcia and Solomon [1983], the 
TEM residual circulation (v*,w*) is computed from the meridional streamfunction (y*). 
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An elliptic equation for x* is obtained from the zonal mean momentum, thermodynamic, 
thermal wind, and continuity equations of the TEM system [e.g., Andrews et ai, 1987], 
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The coefficients (C ZZ ,C Z , C zy , C y , C yy , and Co) depend on empirically determined u 
and T fields (discussed below) and are similar to those derived in Garcia and Solomon 
[1983]. The forcing term (Cf) is proportional to the vertical gradient of the total 
momentum forcing and the latitudinal gradient of the total heating rate, 


Cf — fi ( Fpw + Few + F^w) + Q L A Qgw + Qpw) (A2) 

Here, mechanical forcing is taken as the E-P flux divergence due to dissipating planetary 
waves ( Fpw ), gravity waves (Few), and equatorial Kelvin waves {Fkw)- Thermal 
forcing is due to the combined effects of the net diabatic heating ( Q D ), latent heating 
( Ql ), and the net eddy heat flux associated with planetary waves ( Qpw ) and gravity 
wave breaking in the mesosphere {Qgw)- These quantities will be described in the 
following sections. 

This methodology allows us allows us to compute the planetary wave momentum 
forcing {Fpw in equation Al) and self-consistent K yy fields directly from 3-D 
meterological analyses, instead of inferred from the residual in the momentum equation. 
This approach differs from previous diagnostic calculations which compute the residual 
circulation from only the thermodynamic and continuity equations [e.g., Shine, 1989; 
Eluszkiewicz et al, 1996]. 

The u and T fields used in the coefficients in equation (Al) were compiled primarily 
from the 17-year (1979-1995) average of NCEP data for 1000-1 mbar, and the CIRA-86 
empirical reference atmosphere [ Fleming et al., 1990] for the mesosphere above 1 
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mbar. In most regions outside the tropics, zonal mean winds were derived from the 
temperatures [ Fleming et al . , 1990]. At low latitudes (±20°) where the thermal wind 
derivation is problematic, and at all latitudes above 80 km where the temperature data 
contained in the CIRA model is sparse, we used climatological zonal wind measurements 
from the high resolution Doppler imager (HRDI) onboard UARS. Below 10 km where 
HRDI does not make measurements, zonal wind at 20°S-20°N was derived from a second 
derivative thermal wind calculation [ Fleming et al, 1990]. The temperature and zonal 
wind fields were adjusted where necessary to maintain exact thermal wind balance in all 
regions, as necessitated by the streamfunction formulation. Figure A.la-b shows these 
fields for January and indicates the well known zonal mean structure during solstice 
conditions, as has been discussed frequently in previous climatological compilations 
[e.g., Randel , 1992; Fleming et al. , 1990]. 

To obtain x* » equation (Al) is solved from pole to pole in finite difference form by 
either successive over-relaxation or the method of Lindzen and Kuo [1969]. The altitude 
range extends from the ground to ~116 km. The model chemistry extends up to 90 km, 
and the 90-116 km layer has only an abbreviated representation of the dynamics and 
acts as a sponge layer for the model transport. For boundary conditions, z)*=0 at the 
poles and u)*=0 at the ground (1013 mbar) and the top boundary (0.00008 mbar) so 
that there is no flow into or out of the model domain (^*=0 at the four boundaries). 
The streamfunction formulation satisfies the condition, /+ 9 9 0 ° o ° pw* cos <f>d(j) - 0. 

To obtain reasonable long lived tracer simulations in the stratosphere, it is necessary 
to impose an additional boundary condition on x* one level above the ground (762 
mb - denoted as X 2 *)- Following Garcia et al. [1992], we calculate X 2 as, 

cos f°° /ri 

X 2 - ~ "jp 2 J p(Fpw + F G w + F KW )dz (A3) 

This expression is derived from the linearized steady state momentum equation, 



29 


— fv* — Fpw + Few + Fkw i and follows from the downward control principle under 
quasi-geostrophic scaling [e.g., Haynes et al., 1991; Rosenlof and Holton, 1993]. 
Equation (A3) is applied poleward of 15°N (°S) with values in the tropics interpolated 
between ±15°. The vertical velocity at 762 mbar obtained from x 2 * depends solely on 
the momentum forcing of the free atmosphere above this pressure surface. Upward 
motion out of the boundary layer generally occurs equatorward of ±50° with downward 
motion into the boundary layer at higher latitudes. To obtain good model-measurement 
comparisons of total ozone and long lived tracers, we weakened the circulation by scaling 
the expression in (A3) by a constant factor of .75 for all latitudes and seasons. We note 
that as an independent check of this method, a very similar circulation was obtained by 
computing \ 2 * via the definition, [e.g., Andrews et al., 1987], 
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using climatological averages of the NCEP reanalyses [ Kalnay et al, 1996]. 

Figure A.lc-d shows w* and v* for January. These fields are generally similar to 
previous model calculations [e.g., Garcia et al, 1992]. Upward motion is seen throughout 
the tropical stratosphere, and in the extratropical middle to upper stratosphere and 
mesosphere in summer. Maximum ascent of 3.6 cm/s occurs in the polar upper 
mesosphere. The winter extratropics are characterized by downward motion throughout 
the troposphere and middle atmosphere. The strong upward motion of the Hadley cell 
is also evident in the tropical troposphere. The meridional velocities are generally weak 
in the stratosphere, but reach maximum values of 6-7 m/s in the upper mesosphere 
reflecting the summer to winter meridional flow pattern induced by gravity wave drag 
(section A.3.b). 

Figure A. 2 shows a heuristic view of the model residual circulation, depicted 
by streamlines, and the simulated stratospheric CH4 distribution shown in color, for 
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January and September. This illustrates the controlling influence of the circulation on 
long lived tracers such as CH 4 . Upward motion is present in the tropical stratosphere 
throughout the year, advecting CH 4 -rich air up from the troposphere. Strong descent 
is evident in the winter polar region of both hemispheres, which brings CH 4 -poor air 
down from the mesosphere. The summer lower stratosphere in both hemispheres is 
characterized by very weak vertical motion. 

A. 2. Heating Rates 

The latent heating rates (Ql) are based on the empirical climatology compiled 
by Newell et ai [19/4], shown in Figure A. 3 a. The diabatic heating rates (Qo) are 
computed following Rosenfield et ai [1994], using climatological temperatures (section 
A.l), ozone, and water vapor. The ozone climatology was constructed as follows. Values 
in the troposphere are based on SBUV data as described in McPeters et ai [1984]. 
SAGE II data (both sunrise and sunset) for 1984-1989 are used in the lower stratosphere 
(100-10 mbar). Nimbus-/ SBUV data for 1979-1989 are used for the upper stratosphere 
(10-0.7 mbar), and SME IR and UV measurements for 1982-1986 are used for the 
mesosphere above 0.7 mbar. The water vapor climatology is based on Oort [1983] for 
the ground to 375 mbar, SAGE II data for the upper troposphere and stratosphere, and 
the reference profile of Remsberg et ai [1989] for the mesosphere above 1 mbar. 

The net diabatic heating rates for January are shown in Figure A. 3 b. These 
values are generally consistent with previous heating rate calculations that utilize 
climatological data [e.g., Shine , 1989; Huang and Smith , 1991; Newman and Rosenfield , 
1997]. Net diabatic heating (cooling) is seen in regions where the observed temperature 
is colder (warmer) than the radiative equilibrium temperature. Largest cooling occurs in 
the winter extratropics due to planetary and gravity wave effects, with the latter being 
dominant in the mesosphere. Gravity wave breaking drives the summer mesosphere to 
be far colder than the radiative equilibrium temperature, leading to large net diabatic 
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heating there. Small net diabatic cooling occurs throughout the troposphere during all 
seasons. Largest values are ~-1.5 K/day which occur in the tropical mid-troposphere. 

Gravity wave breaking in the mesosphere induces a net thermal flux due to the 
combination of an eddy heat flux convergence and the turbulent diffusion of the mean 
thermal field. We account for the net heating due to these processes following Schoeberl 
et al. [1983] and Huang and Smith [1991]. Here we utilize the observed temperature 
climatology (section A.l) and vertical diffusivity computed from the gravity wave 
parameterization (section A.3.b). Our calculated net heating from these processes 
{Qgw, Figure A. 3c) is generally similar to that obtained in the previous studies. Net 
cooling occurs throughout almost the entire mesosphere, with a maximum of -12.5 
K/day at 80-90 km in the mid-latitude SH during January. 

Planetary waves also induce a net thermal flux ( Qpw in equation A2). This quantity 
(not shown) can be computed from 3-D meteorological fields [e.g., Andrews et al, 1987], 
but is small compared to the other heating rate contributions. Maximum values are 
generally less than ±.75 K/day occurring in the high latitude winter stratosphere-lower 
mesosphere, and the mid-high latitude troposphere throughout the year. 

The total heating rate for January is shown in Figure A. 3d. Diabatic heating/cooling 
dominates the total throughout the stratosphere, with the gravity wave effects playing 
a significant role in the mesosphere along with the diabatic contribution. Both the 
diabatic and latent heating are important in the troposphere, and the combination of 
these processes leads to horizontal gradients in the total heating rate which determines 
the classical Hadley cell circulation in the tropical troposphere and lower stratosphere 
(Figure A.lc). 

A. 3. Wave Parameterizations 

A. 3. a. Planetary and Synoptic Scale Waves 

Planetary wave dissipation generates rapid irreversible mixing of tracers which 
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can be expressed as a horizontal diffusive transport in two-dimensional models. This 
wave dissipation also induces a drag on the zonal mean flow which can modify the 
mean advective transport by the meridional circulation. Previous investigations have 
illustrated the importance of accounting for these processes in a self-consistent manner 
[Newman et al., 1988; Garcia , 1991]. 

The E-P flux divergence computed from 3-dimensional meteorological analyses 
[e.g., Andrews et al., 1987] provides a diagnostic estimate of the planetary wave drag 
( Fpw in equation A2). Following Randel and Garcia [1994], horizontal mixing rates can 
be calculated self-consistently as the ratio of the E-P flux divergence to the latitudinal 
gradient of zonal mean potential vorticity, K yy = —Fpw/q y . This formulation follows 
from the flux-gradient relationship [e.g., Newman et al., 1988], v T q r = -K yy q y , where the 
E-P flux divergence is taken as the zonal mean PV flux for quasi- geostrophic planetary 
waves. 

For the troposphere and stratosphere, we have derived climatological E-P flux 
divergence and K yy quantities using 17 years (1979-1995) of daily 3-D NCEP analyses 
for 1000-1 mb. This computation includes the wave dissipation and induced mixing 
due to both thermal damping and wave breaking. To account for the effects from both 
planetary and synoptic scale waves in the troposphere and very lower stratosphere, we 
include zonal wavenumbers 1-12 from the ground to 10 mb. Above this level, zonal 
waves 1-6 are included since only the planetary scale waves are important in this region. 

For the eddy winds in the E-P flux calculations, we used balanced winds computed 
from daily NCEP geopotential height analyses for regions poleward of ±20° [Randel, 
1992], and HRDI winds for the tropics between ±20°. Because of data gaps and the 
limited longitudinal coverage of the daily HRDI data, we used monthly mean HRDI 
winds and NCEP temperatures for 1992 to 1996 to compute the tropical E-P flux and 
K yy quantities. Using monthly mean data in this way accounts for the stationary wave 
components but not the transient eddy contributions. However, the HRDI-based wave 
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drive and diffusion quantities resulted in model simulated tracer fields in the tropics and 
subtropics that were in good agreement with observations [ Fleming et al., 1999]. 

For the mesosphere above 1 mbar, Fpw and K yy were computed from the CIRA-86 
stationary planetary wave climatology [ Barnett and Labitzke , 1990]. This contains 
monthly mean temperature and geopotential height data for zonal waves 1-2 based on 
Nimbus-6 satellite measurements for 1975-1978. As in the stratosphere, balanced winds 
were used for the eddy winds poleward of ±20°, with a second derivative balanced wind 
calculation used in the tropics [Fleming et al., 1990]. Unlike the stratosphere, we did 
not utilize the HRDI winds in the tropical mesosphere since using different time periods 
for the eddy temperatures (1975-1978) and eddy winds (1992-1996) was thought to be 
unsatisfactory for these calculations. 

Consistent with previous investigations [Newman et al., 1988], our computations 
revealed mostly negative areas of E-P flux divergence, but also some positive values, 
especially at high latitudes. This implies a wave-induced acceleration of the mean flow, 
and a negative mixing rate via the flux-gradient relationship. While regions of positive 
Fpw may be caused by limitations of the original data sets, some areas of positive E-P 
flux may be real, and the implied negative K yy indicates a breakdown in the assumptions 
underlying the flux-gradient relationship. This is also indicative of the ambiguity in 
interpreting E-P diagnostics for wave propagation and dissipation [e.g., Andrews , 1987; 
Newman et al., 1988]. 

To avoid computing negative mixing rates in practice, regions of positive E-P 
flux divergence were set to zero, implying a zero K yy . The resultant F PW fields were 
smoothed once in the latitude-height domain before being applied in the calculations 
for x and K yy . To ensure that K yy remains well behaved, we also apply the criterion 
that q y > 0.5x10 's 1 as suggested by Randel and Garcia [1994]. This also avoids 
negative mixing rates and implies that large positive values of K yy are caused by eddy 
processes via the E-P flux divergence and not by instability of the zonal mean flow as 
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represented by small or negative q y . 

We impose a very small background K yy of 10 8 cm 7 s~ l throughout the stratosphere 
and mesosphere. In the low-middle troposphere, we set a large background K yy of 
2xl0 10 cm 2 s~ 1 to simulate the interhemispheric transit time of 1-2 years estimated from 
SF 6 and halocarbon data near the ground [e.g., Geller et al, 1997], We also impose 
a small K yy of 10 8 cm 2 s 1 around the tropopause region where the tropopause height 
changes across adjacent latitudes. This decreased the anomolous cross-tropopause 
transport and gave better model simulations of water vapor compared with observations 
in the very lower stratosphere at middle to high latitudes. Because of the uncertainties 
involved and the relatively small effect, we do not include the transport due to the 
off-diagonal K yz and K zy terms in the diffusion tensor. These arise from the projection 
of K yy from isentropes to pressure surfaces [Newman et al., 1988]. 

Figure A. 4 shows the planetary wave driving ( Fpw ) for January. In the troposphere, 
large wave drive is evident poleward of ±20° during winter and to a lesser extent, 
summer, due to synoptic and planetary wave activity. These features are similar in 
both hemispheres (July not shown). Large Fpw also occurs in the winter extra-tropical 
stratosphere and mesosphere. 

The corresponding stratospheric K yy fields for January and September are shown 
in color in Figure A. 5, and illustrate the regions of large mixing associated with large 
planetary wave activity in the winter middle atmosphere. Here we have also overlaid the 
model CH4 values to show the effect of strong mixing on the long lived tracer fields, for 
example as seen in the winter midlatitude surf-zone region in the stratosphere. We note 
that consistent with Randel and Garcia [1994], the direct ratio calculation of K yy gives 
large mixing rates at high latitudes in the upper stratosphere and lower mesophere, and 
also at 20°-45°N above 100 mb coincident with weak q y in the surf-zone region. This 
is in contrast to Newman et al. [1988] who computed largest mixing rates only at high 
latitudes using a linear regression fit of Vq' and q y with a zero intercept. Weak mixing 
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is evident in the tropics and summer hemisphere, except for the very lower summer 
stratosphere where decaying baroclinic wave activity from the troposphere generates 
large mixing below 20 km. 

The seasonal variation of K yy and u at several different levels are shown in Figure 
A. 6. In the mid-upper stratosphere and mesosphere, weak mixing occurs in the easterly 
flow regimes of the summer hemisphere and in the tropics, consistent with a minimum 
of planetary wave activity in these regions. Note that the tropical stratospheric 
mixing rates based on the HRDI winds are weak but not negligible, with values as 
large as 8xl0 9 cm 2 s~ 1 . In the mid-late winter, large K yy values are seen along the 
subtropical edges of the polar night jets, characterizing the surf-zone region of the winter 
stratosphere. This area of maxima moves poleward during late winter, with large mixing 
at high latitudes during spring of both hemispheres associated with the deceleration and 
break up of the polar vorticies. These figures also show the contrasting degree of winter 
polar vortex isolation in the two hemispheres, reflecting the hemispheric asymmetry in 
planetary wave activity. Weak mixing occurs in the core of the SH mid-winter vortex, 
consistent with the inhibition of planetary wave propagation through strong westerly 
winds [Charney and Drazin , 1961]. The NH exhibits significantly weaker u and larger 
Kyy throughout the winter stratosphere poleward of 25°N. In the lower mesosphere, both 
hemispheres show strong subtropical jets near ±35° in mid-winter, with strong mixing 
occurring mainly on the poleward flanks of the jets. In the very lower stratosphere 
(68 mb), westerly winds and significant mixing occur throughout much of the year in 
the extratropics of both hemispheres. There is also evidence of a semi-annual signal 
in the mixing in both hemispheres, with largest K yy in spring and fall coinciding with 
the build-up and decrease of the winter polar jets. As at higher levels, weaker mixing 
at 68 mbar is observed in the tropics, and in mid-winter in the cores of the polar jets, 
especially in the SH. 
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A.3.b. Gravity Waves 

Computations of vertical eddy diffusion ( K zz ) and mechanical forcing {Few in 
equations A1-A2) from gravity waves are based on the parameterization originally 
formulated by Lindzen [1981] and modified by Holton and Zhu [1984] (hereafter HZ84). 
Other studies have used this parameterization to compute the time evolution of the 
zonal mean flow [ Schoeberl et al., 1983; HZ84; Garcia and Solomon, 1985], However, 
since u and T are pre-determined in our model, we compute the seasonal and spatial 
distribution of wave drag and diffusion based on the empirically determined zonal 
wind and temperature fields and an input set of gravity wave parameters. These 
computed wave effects are highly dependent on the background wind field in the 
middle atmosphere [e.g., Shine , 1989], To account for the longitudinal variations of the 
background wind which can be significant during the NH winter, we use a 3-D zonal 
wind field in the parameterization. This wind field is composed of the zonal mean 
wind described in section A.l, with the longitudinal variability based on NCEP and 
HRDI data as discussed in section A. 3. a. The longitudinal grid resolution is 30° which, 
although rather coarse, should be sufficient to accomodate the planetary scale features 
of the middle atmosphere. 

In accordance with the tropospheric sources of gravity waves, we assume a spectrum 
of 9 waves with zonal phase speeds c = 0, ±10, ±20, ±30, ±40 ms~ x , all with a horizontal 
wavelength of 200 km (HZ84). Here we have assumed that waves are only traveling 
zonally, although similar results are obtained by assuming an isotropic spectrum of 
waves as in HZ84. Following Garcia et al. [1992], we include two additional waves in the 
tropics with phase speeds c = ±50 ms -1 and zonal wavelength of 800 km. At the lower 
boundary (375 mbar), we specify the vertical velocity (u3o) as a function of the wave 
phase speed, with the slowest phase velocity waves having the largest w 0 . The velocities 
range from .03 ms 1 for the slowest waves to .002 ms -1 for the fastest waves, with these 
values skewed to have slightly larger values for the westerly waves [Huang and Smith, 
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1991]. We impose only a weak geographical variation on w 0 , with slightly larger values 
at midlatitudes compared with the tropics and polar regions. To mimic the generation 
of gravity waves by convection, we also specify slightly larger w 0 over subtropical land 
areas during summer following Wu and Waters [1996]. However, w r e note that in our 
current model formulation, the geographical and seasonal variability in the computed 
gravity wave effects is determined primarily by the background wind. 

For each latitude and longitude, the algorithm determines the breaking level z^, 
turbulent diffusion, and vertical momentum flux convergence for each w r ave component. 
The contributions to diffusion and wave drag are summed over all wave components to 
obtain distributions of K zz and Fgw ■ We also include in Few the wave driving due 
to the turbulent diffusion of zonal momentum. The K zz and F G w values are zonally 
averaged, and smoothed in latitude and height. Following previous studies, we apply an 
efficiency factor (e = .025) to K zz and F G w to account for the fact that gravity wave 
breaking is not continuous, but rather is sporadic in time and space. The values of e and 
wq (used to determine z^r ) were specified to compute a wave drag peaking at about 100 
ms l day 1 in the upper mesosphere, in rough qualitative agreement with that inferred 
previously [e.g., Holton, 1983; Fritts and Vincent, 1987; Shine, 1989; Huang and Smith, 
1991]. Our computed K zz and Few gave good agreement in the computed seasonal cycle 
of mesospheric H 2 O compared with UARS/HALOE measurements (see Part 1). It may 
be possible to improve the model simulations by tuning the seasonal and geographical 
variability of e and xvq as observational characterization of gravity waves improve. 

We also include a contribution to the zonal mean F cw and K zz from the breaking 
diurnal tide in the lower thermosphere. In accordance with UARS/HRDI observations 
[Lieberman and Hays, 1994] and theoretical estimates [Forbes, 1982], w'e specify the 
largest values at the equator at 90-100 km exponentially decreasing in latitude and 
height, w-ith a semiannual variation maximizing at the equinoxes. Largest values for 
March at 90 km at the equator are set to -10 ms~^ day~ l for wave drag and 30 m 2 s~* 
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for diffusion. 

The total zonal mean F G w from gravity waves and the diurnal tide is then applied 
in equations A1-A2. A background K zz , specified to increase from its value at 30 km 
(section A.3.c) to .3 m 2 s _1 in the mesosphere, is added to that due to gravity waves 
and the diurnal tide. The total K zz is then used in the model transport of chemical 
constituents. We impose an upper limit of 100 m 2 s _1 to accomodate the 2.4 hour time 
step used for the vertical diffusive transport of constituents. 

Figure A. 7 shows latitude-height sections of wave drag and vertical diffusion 
from the gravity wave scheme, diurnal tidal effects, prescribed background, and 
troposphere/lower stratosphere specification (see section A.3.c). These calculations are 
consistent with previous model results [e.g., Garcia and Solomon , 1985], The largest 
wave drag occurs in the middle to upper mesosphere. Values of F G w are opposite in 
sign to u (Figure A. lb) as the gravity wave drag acts to close off the jets in the upper 
mesosphere, and induce the summer- winter meridional flow pattern (Figure A. Id). The 
K zz field also shows maxima occurring in regions of large u in the mid-latitude upper 
mesosphere, with a distinct minimum in the tropics. This pattern is very similar to 
zonally averaged small scale variances seen in UARS microwave limb sounder (MLS) 
radiance measurements, which may be an indicator of gravity wave activity [ Wu and 
Waters, 1996]. The tropopause is also evident in Figure A.7b, as seen by the sharp 
transition from large K zz in the troposphere to very small values in the stratosphere. 

A.S.c. Tropospheric/Stratospheric K zz 

In the troposphere and lower stratosphere below 30 km, K zz is based on the 
climatological temperature lapse rate, and is determined as follows. We specify the 
largest K zz to be 100 m 2 s _1 for the steepest lapse rate observed in the lower troposphere 
(-8 K/km). Following the analyses of Mote et al. [1998] and Hall and Waugh [1997], we 
specify the minimum K zz to be .01 m 2 s _1 in the tropical lower stratosphere where the 
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lapse rate is most stable (~ +AK/km). The K zz values at each grid point throughout 
the year are then interpolated from these extreme values using the corresponding lapse 
rate. This methodology assumes that larger mixing occurs in regions of weaker static 
stability, given that the latter is indicative of stronger convective overturning. Here 
we have not accounted for the effects induced by instability in the vertical wind shear 
given the difficulty in quantifying the diffusion due to this process. As seen in Figure 
A. 7b, this K zz profile mimics the large vertical mixing in the troposphere with a sharp 
decrease in mixing across the tropopause into the lower stratosphere. 

A.S.d. Kelvin Waves 

Equatorial Kelvin waves are thought to be important in driving the westerly 
acceleration of the stratospheric SAO and related meridional circulation. To approximate 
this process, we used the parameterization of Dunkerton [1979] and Gray and Pyle 
[1987] along with the empirical u and T (section A.l) to diagnose the mechanical forcing 
from thermally damped Kelvin waves ( F KW in equations A1-A2). We assume a zonal 
wave phase speed of 50 ms~ i for zonal wavenumber 1. The vertical momentum flux 
(B) at the bottom boundary (138 mbar) at the equator is specified to be .0105 m 2 s~ 2 
which is a factor of 1.5 larger than that used by Gray and Pyle [1987]. We also used a 
thermal damping rate (a) which was slightly slower than the ‘fast’ rates of Dunkerton 
[1979], with a minimum e-folding damping time of 2.7 days near 50 km. We found that 
specifying these values of B and a resulted in model simulated tracer fields most similar 
to observations. 

Figure A. 8 show u and the computed F^w quantity at the equator. The zonal 
mean wind shows the well known SAO signature above 10 mb, with a phase reversal 
between the upper stratosphere and upper mesosphere. Maximum F K w occurs near 
the stratopause, with values of 2-3 ms~ l day~ 1 . There is a semiannual oscillation in 
Fkwi with significant values extending into the upper mesosphere during January 



and July, coincident with the weakest westerly u in the lower mesosphere. Weakest 
westerly accelerations in the upper stratosphere and lower mesosphere occur during the 
equinoxes, when the westerlies are at a maximum. This is consistent with the modeling 
study of Dunkerton [1979] and the observations of Hirota [1978]. 

Finally, we show the total wave driving ( Ftot = Fpw + Few + Fkw) for 
January in Figure A. 9. Gravity w r ave drag dominates the total throughout the upper 
mesosphere, with the planetary wave effects being a greater contribution to the total in 
the stratosphere. The signal in the middle to high latitude troposphere is due to the 
combination of planetary and baroclinic eddies. 
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Figure Captions 


Figure 1. Vertical profiles of carbon 14 at 31°N for 
July 1964 and July 1966. Plotted are the time depen- 
dent model simulations for base scenario A (solid line) 
and the four circulation scenarios B-E listed in Table 
1: very weak circulation (B, dashed-triple dot line); 
weak circulation (C, dashed line); strong circulation (D, 
dashed-dot line); and very strong circulation (E, dotted 
line). Also shown are the observations from Kinnison 
et al. [1994] (dashed-asterisk line). Values are in mixing 
ratio units, defined as 10 5 atoms of 14 C per gram of dry 
air. 


Figure 1. Vertical profiles of carbon 14 at 31°N for July 1964 and July 1966. Plotted are 
the time dependent model simulations for base scenario A (solid line) and the four circulation 
scenarios B-E listed in Table 1: very weak circulation (B, dashed-triple dot line); weak circulation 
(C, dashed line); strong circulation (D, dashed-dot line); and very strong circulation (E, dotted 
line). Also shown are the observations from Kinnison et al. [1994] (dashed-asterisk line). Values 
are in mixing ratio units, defined as 10 5 atoms of 14 C per gram of dry air. 

Figure 2. Annually averaged mean age of air (years) 
derived from time dependent model simulations of SF6 
for the base scenario A and the four circulation sce- 
narios B-E listed in Table 1. The age is taken relative 
to the global mean value at the surface. The contour 
interval is 1 year. 

Figure 2. Annually averaged mean age of air (years) derived from time dependent model 
simulations of SF 6 for the base scenario A and the four circulation scenarios B-E listed in Table 
1. The age is taken relative to the global mean value at the surface. The contour interval is 1 
year. 

Figure 3. Age of air derived from SF& as a function of 
latitude for October/November 1994 using ER-2 mea- 
surements at 19-21 km (triangles) and the model cir- 
culation scenarios listed in Table 1 at 20 km (lines). 

Included are the base model scenario (A, solid line), 
the very weak circulation (B, dashed-triple dot line), 
the weak circulation (C, dashed line), the strong circu- 
lation (D, dashed-dot line), and the very strong circu- 
lation (E, dotted line). The age is taken relative to the 
global mean value at the surface. 

Figure 3. Age of air derived from SF 6 as a function of latitude for October/November 1994 
using ER-2 measurements at 19-21 km (triangles) and the model circulation scenarios listed in 
Table 1 at 20 km (lines). Included are the base model scenario (A, solid line), the very weak 
circulation (B, dashed-triple dot line), the weak circulation (C, dashed line), the strong circulation 
(D, dashed-dot line), and the very strong circulation (E, dotted line). The age is taken relative 
to the global mean value at the surface. 



Figure 4. Vertical profiles of age of air derived from 
SFg from the model circulation scenarios listed in Table 
1 and balloon data taken at the latitudes and seasons 
indicated. The observations are from OMS data (trian- 
gles) and the balloon data of Hamisch et al [1996] (as- 
terisks). Included are the base model scenario (A, solid 
line), the very weak circulation (B, dashed-triple dot 
line), the weak circulation (C, dashed line), the strong 
circulation (D, dashed-dot line), and the very strong 
circulation (E, dotted line). The age is taken relative 
to the global mean value at the surface. 


Figure 4. Vertical profiles of age of air derived from SF 6 from the model circulation scenarios 
listed in Table 1 and balloon data taken at the latitudes and seasons indicated. The observations 
are from OMS data (triangles) and the balloon data of Hamisch et al. [1996] (asterisks). Included 
are the base model scenario (A, solid line), the very weak circulation (B, dashed-triple dot line), 
the weak circulation (C, dashed line), the strong circulation (D, dashed-dot line), and the very 
strong circulation (E, dotted line). The age is taken relative to the global mean value at the 
surface. 


Figure 5. Vertical profiles of CO 2 (ppmv) versus N 2 O 
(ppbv) for July/August 1994 at 70°S to 60°S from the 
ER-2 measurements (crosses) and the model circulation 
scenarios listed in Table 1. Included are the base model 
scenario (A, solid line), the very weak circulation (B, 
dashed-triple dot line), the weak circulation (C, dashed 
line), the strong circulation (D, dashed-dot line), and 
the very strong circulation (E, dotted line). Approxi- 
mate altitude is shown on the right-hand axes. 


Figure 5. Vertical profiles of CO 2 (ppmv) versus N 2 O (ppbv) for July/August 1994 at 70°S to 
60°S from the ER-2 measurements (crosses) and the model circulation scenarios listed in Table 1. 
Included are the base model scenario (A, solid line), the very weak circulation (B, dashed-triple 
dot line), the weak circulation (C, dashed line), the strong circulation (D, dashed-dot line), and 
the very strong circulation (E, dotted line). Approximate altitude is showm on the right-hand 
axes. 


Figure 6. Vertical profiles of carbon 14 at 31°N for 
July 1966. Plotted are the time dependent model sim- 
ulations for base scenario A (solid line) and the three 
scenarios (F-H) of stratospheric minimum K zz listed in 
Table 1: minimum A\ 2 =0.001 m 2 /s (F, dashed-triple 
dot line); minimum K zz ~ 0.1 m 2 /s (G, dashed line); and 
minimum m 2 /s (H, dotted line). Also shown are 

the observations from Kinnison et al [1994] (dashed- 
asterisk line). Values are in mixing ratio units, defined 
as 10 5 atoms of 14 C per gram of dry air. 


Figure 6. Vertical profiles of carbon 14 at 31°N for July 1966. Plotted are the time dependent 
model simulations for base scenario A (solid line) and the three scenarios (F-H) of stratospheric 
minimum K zz listed in Table 1: minimum K zz ~0.001 m 2 /s (F, dashed-triple dot line); minimum 
K zz = 0.1 m 2 /s (G, dashed line); and minimum K zz = 1 m 2 /s (H, dotted line). Also shown are the 
observations from Kinnison et al. [1994] (dashed-asterisk line). Values are in mixing ratio units, 
defined as 10 5 atoms of 14 C per gram of dry air. 



Figure 7. Annually averaged mean age of air (years) 
derived from time dependent model simulations of Sp6 
for the base scenario A and the three scenarios (F-H) of 
stratospheric minimum K zz listed in Table 1. The age 
is taken relative to the global mean value at the surface. 
The contour interval is 1 year. 


Figure 7. Annually averaged mean age of air (years) derived from time dependent model 
simulations of SF6 for the base scenario A and the three scenarios (F-H) of stratospheric minimum 
K zz listed in Table 1. The age is taken relative to the global mean value at the surface. The 
contour interval is 1 year. 


Figure 8. Age of air derived from SFg as a function of 
latitude for October/November 1994 using ER-2 mea- 
surements at 19-21 km (triangles) and the model sce- 
narios of stratospheric minimum K zz listed in Table 1 
at 20 km (lines). Included are the base model scenario 
(A, solid line), minimum K zz =0.001 m 2 /s (F, dashed- 
triple dot line); minimum K zz — 0.1 m 2 /s (G, dashed 
line); and minimum K zz — 1 m 2 /s (H, dotted line). The 
age is taken relative to the global mean value at the 
surface. 


Figure 8. Age of air derived from SF6 as a function of latitude for October/November 1994 using 
ER-2 measurements at 19-21 km (triangles) and the model scenarios of stratospheric minimum 
K zz listed in Table 1 at 20 km (lines). Included are the base model scenario (A, solid line), 
minimum #--=0.001 m 2 /s (F, dashed-triple dot line); minimum #**=0.1 m 2 /s (G, dashed line); 
and minimum #** = 1 m 2 /s (H, dotted line). The age is taken relative to the global mean value 
at the surface. 

Figure 9. Vertical profiles of age of air derived from 
SF 6 from the model scenarios of stratospheric minimum 
K zz listed in Table 1 and balloon data taken during 
September near 40°N. The observations are from OMS 
data (triangles) and the balloon data of Hamisch et 
al. [1996] (asterisks). Included are the base model sce- 
nario (A, solid line), the minimum #**=0.001 m 2 /s (F, 
dashed-triple dot line); minimum K zz — 0.1 m 2 /s (G, 
dashed line); and minimum #** = 1 m 2 /s (H, dotted 
line). The age is taken relative to the global mean value 
at the surface. 

Figure 9. Vertical profiles of age of air derived from SF6 from the model scenarios of strato- 
spheric minimum K zz listed in Table 1 and balloon data taken during September near 40°N. 
The observations are from OMS data (triangles) and the balloon data of Hamisch et al [1996] 
(asterisks). Included are the base model scenario (A, solid line), the minimum K zz =0.001 m 2 /s 
(F, dashed-triple dot line); minimum #\*=0.1 m 2 /s (G, dashed line); and minimum K zz — 1 m 2 /s 
(H, dotted line). The age is taken relative to the global mean value at the surface. 



Figure 10. Seasonal profiles of zonal and monthly 
mean H 2 O from the HALOE version 19 data (dashed- 
asterisk) and three model scenarios at 80 km (.01 mbar) 
for 45°N, equator, and 45°S. Included are the base sce- 
nario (A, solid line), and scenarios with weak (I, dashed 
line) and strong (J, dotted line) mesospheric gravity 
wave effects. 


Figure 10. Seasonal profiles of zonal and monthly mean H 2 O from the HALOE version 19 data 
(dashed-asterisk) and three model scenarios at 80 km (.01 mbar) for 45°N, equator, and 45°S. 
Included are the base scenario (A, solid line), and scenarios with weak (I, dashed line) and strong 
(J, dotted line) mesospheric gravity wave effects. 


Figure 11 . Annually averaged mean age of air (years) 
derived from time dependent model simulations of SFq 
for the base scenario A and the 7 K yy scenarios (K- 
Q) listed in Table 1. In panels a-d and g-h, the K yy 
and wave drive are changed above the tropopause. Sce- 
nario (Q) in panel f is the same as the base case (sce- 
nario A, panel e), except K yy has been set to a very 
small constant value of 10 8 cm 2 /s throughout the tropi- 
cal stratosphere (15°S-15°N, 15-50 km) for all seasons. 
In all panels, the age is taken relative to the global mean 
value at the surface. The contour interval is 1 year. 


Figure 11. Annually averaged mean age of air (years) derived from time dependent model 
simulations of SF6 for the base scenario A and the 7 K yy scenarios (K-Q) listed in Table 1. In 
panels a-d and g-h, the K yy and wave drive are changed above the tropopause. Scenario (Q) in 
panel f is the same as the base case (scenario A, panel e), except K yy has been set to a very 
small constant value of 10 8 cm 2 /s throughout the tropical stratosphere (15°S-15°N, 15-50 km) for 
all seasons. In all panels, the age is taken relative to the global mean value at the surface. The 
contour interval is 1 year. 


Figure 12. Age of air derived from SF 6 as a func- 
tion of latitude for October/November 1994 using ER- 
2 measurements at 19-21 km (triangles) and several 
model Kyy scenarios listed in Table 1 at 20 km (lines). 
Included are the base model scenario (A, solid line), 
K yy reduced by a factor of 5 above the tropopause (K, 
dashed-triple dot line), K vy reduced by a factor of 2 
above the tropopause (L, dashed line), K yy increased by 
a factor of 2 above the tropopause (M, dashed-dot line), 
K yy increased by a factor of 5 above the tropopause (N, 
dotted line), and the base model with K yy set to a very 
small constant value of 10 8 cm 2 /s throughout the trop- 
ical stratosphere (15°S-15°N, 15-50 km) for all seasons 
(Q, sold-asterisk line). The age is taken relative to the 
global mean value at the surface. 


Figure 12. Age of air derived from SF 6 as a function of latitude for October/November 1994 
using ER-2 measurements at 19-21 km (triangles) and several model K yy scenarios listed in Table 
1 at 20 km (lines). Included are the base model scenario (A, solid line), K yy reduced by a factor 
of 5 above the tropopause (K, dashed-triple dot line), K yy reduced by a factor of 2 above the 
tropopause (L, dashed line), K yy increased by a factor of 2 above the tropopause (M, dashed-dot 
line), K yy increased by a factor of 5 above the tropopause (N, dotted line), and the base model 
with K vy set to a very small constant value of 10 8 cm 2 /s throughout the tropical stratosphere 
(15°S-15°N, 15-50 km) for all seasons (Q, sold-asterisk line). The age is taken relative to the 
global mean value at the surface. 



Figure 13. Vertical profiles of age of air derived from 
SF 6 from several model K yy scenarios listed in Table 
1, and balloon data taken at the latitudes and seasons 
indicated. The observations are from OMS data (trian- 
gles) and the balloon data of Harnisch et al [1996] (as- 
terisks). Included are the base model scenario (A, solid 
line), K yy reduced by a factor of 5 above the tropopause 
(K, dashed-triple dot line), K yy reduced by a factor of 2 
above the tropopause (L, dashed line), K yy increased by 
a factor of 2 above the tropopause (M, dashed-dot line), 
and Kyy increased by a factor of 5 above the tropopause 
(N, dotted line). The age is taken relative to the global 
mean value at the surface. 


Figure 13. Vertical profiles of age of air derived from SF $ from several model K vy scenarios 
listed in Table 1, and balloon data taken at the latitudes and seasons indicated. The observations 
are from OMS data (triangles) and the balloon data of Harnisch et ai [1996] (asterisks). Included 
are the base model scenario (A, solid line), K yy reduced by a factor of 5 above the tropopause 
(K, dashed-triple dot line), K yy reduced by a factor of 2 above the tropopause (L, dashed line), 
Kyy increased by a factor of 2 above the tropopause (M, dashed-dot line), and K yy increased by 
a factor of 5 above the tropopause (N, dotted line). The age is taken relative to the global mean 
value at the surface. 


Figure 14. Age of air derived from SF 6 as a func- 
tion of latitude for October/November 1994 using ER- 
2 measurements at 19-21 km (triangles) and several 
model Kyy scenarios listed in Table 1 at 20 km (lines). 
Included are the base model scenario (A, solid line), 
and the scenarios with K yy reduced by a factor of 5 
above the tropopause (K, dashed-triple dot line) and 
Kyy increased by a factor of 5 above the tropopause 
(N, dotted line) as in Figure 12. Also included are 
scenarios in which both the K yy and wave drive have 
been reduced by a factor of 5 above the tropopause (O, 
dashed-dot line), and increased by a factor of 5 above 
the tropopause (P, dashed line). The age is taken rela- 
tive to the global mean value at the surface. 


Figure 14. Age of air derived from SF 6 as a function of latitude for October/November 1994 
using ER-2 measurements at 19-21 km (triangles) and several model K vy . scenarios listed in Table 
1 at 20 km (lines). Included are the base model scenario (A, solid line), and the scenarios with 
Kyy reduced by a factor of 5 above the tropopause (K, dashed-triple dot line) and K yy increased 
by a factor of 5 above the tropopause (N, dotted line) as in Figure 12. Also included are scenarios 
in which both the K yy and wave drive have been reduced by a factor of 5 above the tropopause 
(O, dashed-dot line), and increased by a factor of 5 above the tropopause (P, dashed line). The 
age is taken relative to the global mean value at the surface. 

Figure 15. Annually averaged mean age of air (years) 
derived from time dependent model simulations of SF 6 . 

Included are the base scenario A and the four combi- 
nation transport scenarios (S-V) listed in Table 1. Also 
included is the simulation from the previous 1995 ver- 
sion of the model transport in panel f (scenario W). In 
all panels, the age is taken relative to the global mean 
value at the surface. The contour interval is 1 year. 

Figure 15. Annually averaged mean age of air (years) derived from time dependent model 
simulations of SFg. Included are the base scenario A and the four combination transport scenarios 
(S-V) listed in Table 1. Also included is the simulation from the previous 1995 version of the 
model transport in panel f (scenario W). In all panels, the age is taken relative to the global 
mean value at the surface. The contour interval is 1 year. 



Figure 16. Age of air derived from SF 6 as a function of 
latitude for October/November 1994 using ER-2 mea- 
surements at 19-21 km (triangles) and several model 
scenarios listed in Table 1 at 20 km (lines). Included 
are the base model scenario (A, solid line), scenario S 
(dashed-triple dot line) which combines scenarios C, F, 
and M to give a moderately long stratospheric residence 
time, scenario T (dotted line) which combines scenar- 
ios D, G, and L to give a moderately short stratospheric 
residence time, and scenario W (dashed line) which uses 
the previous 1995 version of the model transport. The 
age is taken relative to the global mean value at the 
surface. 


Figure 16. Age of air derived from SFg as a function of latitude for October/November 1994 
using ER-2 measurements at 19-21 km (triangles) and several model scenarios listed in Table 1 
at 20 km (lines). Included are the base model scenario (A, solid line), scenario S (dashed-triple 
dot line) which combines scenarios C, F, and M to give a moderately long stratospheric residence 
time, scenario T (dotted line) which combines scenarios D, G, and L to give a moderately short 
stratospheric residence time, and scenario W (dashed line) which uses the previous 1995 version 
of the model transport. The age is taken relative to the global mean value at the surface. 


Figure 17. Vertical profiles of age of air derived from 
SF 6 from several model scenarios listed in Table 1, and 
balloon data taken at the latitudes and seasons indi- 
cated. The observations are from OMS data (triangles) 
and the balloon data of Hamisch et al [1996] (aster- 
isks). Included are the base model scenario (A, solid 
line), scenario S (dashed-triple dot line) which com- 
bines scenarios C, F, and M to give a moderately long 
stratospheric residence time, scenario T (dotted line) 
which combines scenarios D, G, and L to give a mod- 
erately short stratospheric residence time, and scenario 
W (dashed line) which uses the previous 1995 version 
of the model transport. The age is taken relative to the 
global mean value at the surface. 


Figure 17. Vertical profiles of age of air derived from SF6 from several model scenarios listed 
in Table 1, and balloon data taken at the latitudes and seasons indicated. The observations are 
from OMS data (triangles) and the balloon data of Hamisch et al. [1996] (asterisks). Included 
are the base model scenario (A, solid line), scenario S (dashed-triple dot line) which combines 
scenarios C, F, and M to give a moderately long stratospheric residence time, scenario T (dotted 
line) which combines scenarios D, G, and L to give a moderately short stratospheric residence 
time, and scenario W (dashed line) which uses the previous 1995 version of the model transport. 
The age is taken relative to the global mean value at the surface. 

Figure 18. Scatter plots of the annual and globally 
averaged mixing ratio (ppbv) of Cl y (top) and N 2 0 
(bottom), versus mean age (years) at 27 km for the 
23 transport scenarios as listed in Table 1. Letters for 
each scenario are color coded according to the type of 
transport variation imposed, as indicated in the legend. 

Figure 18. Scatter plots of the annual and globally averaged mixing ratio (ppbv) of Cl y (top) 
and N 2 0 (bottom), versus mean age (years) at 27 km for the 23 transport scenarios as listed in 
Table 1. Letters for each scenario are color coded according to the type of transport variation 
imposed, as indicated in the legend. 



Figure 19. Scatter plots of the annual, global, and 
stratospherically averaged (pressure weighted) lifetime 
(years) of CFCI3 (top) and N 2 0 (bottom), versus mean 
age (years) for the 23 transport scenarios listed in Table 
1. Letters for each scenario are color coded according to 
the type of transport variation imposed, as indicated in 
the legend. Also included are the best empirical lifetime 
estimates (solid line) and the correspoding 1 standard 
deviation uncertainties (dashed lines) from Volk et al. 
[1997]. 


Figure 19. Scatter plots of the annual, global, and stratospherically averaged (pressure weighted) 
lifetime (years) of CFCI3 (top) and N 2 0 (bottom), versus mean age (years) for the 23 transport 
scenarios listed in Table 1. Letters for each scenario are color coded according to the type of 
transport variation imposed, as indicated in the legend. Also included are the best empirical 
lifetime estimates (solid line) and the correspoding 1 standard deviation uncertainties (dashed 
lines) from Volk et al [1997]. 

Figure A.l. Latitude-height sections of zonal mean 
temperature, zonal wind, and residual vertical velocity 
and horizontal velocity from the model transport. The 
contour intervals are 10 K for temperature; 10 m/s for 
zonal wind; 0.5 cm/s for w* including the contours for 
±.2,±.l, and ±.05 cm/s; and 2 m/s for v* including 
the contours for ±1, ±.5, and ±.2 m/s. Negative values 
are shaded. 

Figure A.l. Latitude-height sections of zonal mean temperature, zonal wind, and residual 
vertical velocity and horizontal velocity from the model transport. The contour intervals are 10 
K for temperature; 10 m/s for zonal wind; 0.5 cm/s for w* including the contours for ±.2,±.l, 
and ±.05 cm/s; and 2 m/s for v* including the contours for ±1,±.5, and ±.2 m/s. Negative 
values are shaded. 

Figure A. 2 . Model simulation of CH 4 in ppbv (color) 
for January and September, along with streamlines 
depicting the sense of the model residual circulation 
(v*,w m ). w* has been multiplied by a deformation fac- 
tor of 200, approximating the ratio of N/ f. 

Figure A. 2. Model simulation of CH 4 in ppbv (color) for January and September, along with 
streamlines depicting the sense of the model residual circulation (v*,w*). w * has been multiplied 
by a deformation factor of 200, approximating the ratio of N/f. 

Figure A. 3. Latent heating rates (a), net diabatic 
heating (b), net gravity wave heating (c), and the to- 
tal heating rate (d) for January. See text for details. 

Contour intervals are .5 K/day for latent heating, and 
2 K/day including the contours for ±1, and ±.5 K/day 
for panels b, c, and d. Negative values are shaded. 

Figure A. 3. Latent heating rates (a), net diabatic heating (b), net gravity wave heating (c), 
and the total heating rate (d) for January. See text for details. Contour intervals are .5 K/day 
for latent heating, and 2 K/day including the contours for ±1, and ±.5 K/day for panels b, c, 
and d. Negative values are shaded. 
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Figure A. 4. Eliassen-Palm flux divergence due to 
planetary and synoptic scale waves (Fpw) f° r January. 

The contour intervals are 2 m/s/day including ±1 and 
±.5 m/s/day (values < -2 m/s/day are shaded). 

Figure A. 4. Eliassen-Palm flux divergence due to planetary and synoptic scale waves (Fpvv) 
for January. The contour intervals are 2 m/s/day including ±1 and ±.5 m/s/day (values < 
-2 m/s/day are shaded). 

Figure A. 5. Model simulation of CH 4 (ppbv, solid 
contours) as in Figure A. 2 for January and September, 
along with the model K yy fields depicted in color. The 
contour interval for CH 4 is 100 ppbv, and the K yy values 
are in 10 8 cm 2 /s. 

Figure A. 5. Model simulation of CH 4 (ppbv, solid contours) as in Figure A. 2 for January and 
September, along with the model K yy fields depicted in color. The contour interval for CH 4 is 
100 ppbv, and the K yy values are in 10 8 cm 2 /s. 

Figure A. 6. Month-latitude sections of the model K yy 
fields in color, along with u overlaid in solid contours, 
for 68 mbar, 12 mbar, 1.7 mbar, and .1 mbar. The 
contour interval for u is 20 m/s, including ±10 m/s, 
and the K yy values are in 10 8 cm 2 /s. 

Figure A. 6. Month-latitude sections of the model K yy fields in color, along with u overlaid 
in solid contours, for 68 mbar, 12 mbar, 1.7 mbar, and .1 mbar. The contour interval for u is 
20 m/s, including ±10 m/s, and the K yy values are in 10 8 cm 2 /s. 

Figure A. 7. Eliassen-Palm flux divergence due to 
gravity waves (Fqw) (a)> and the associated model K zz 
values (b) for January. The contour intervals are, for 
(a): 20 m/s/day including ±10, ±5, and ±1 m/s/day 
(negative values are shaded); (b): 20 m 2 /s including 
.02, .03, .1, 1, 5, 10 m 2 /s. 

Figure A. 7. Eliassen-Palm flux divergence due to gravity waves (Few) (a), and the associated 
model K zz values (b) for January. The contour intervals are, for (a): 20 m/s/day including 
±10, ±5, and ±1 m/s/day (negative values are shaded); (b): 20 m 2 /s including .02, .03, .1, 1, 5, 
10 m 2 /s. 

Figure A. 8. Month-height sections of u and the 
Eliassen-Palm flux divergence due to Kelvin waves 
(Fkw) at the equator. The contour intervals are 
for 10 m/s for u (negative values are shaded), and 
.5 m/s/day including .1 m/s/day for Fkw • 

Figure A. 8. Month-height sections of u and the Eliassen-Palm flux divergence due to Kelvin 
waves {Fkw) at the equator. The contour intervals are for 10 m/s for u (negative values are 
shaded), and .5 m/s/day including .1 m/s/day for Fkw- 

Figure A. 9. The total Eliassen-Palm flux diver- 
gence ( Ftot ) for January. Contour intervals are 
20 m/s/day and include the contours for ±1,±2,±5, 
and ±10 m/s/day. Negative values are shaded. 

Figure A. 9. The total Eliassen-Palm flux divergence {Ftot) for January. Contour intervals 
are 20 m/s/day and include the contours for ±1, ±2, ±5, and ±10 m/s/day. Negative values are 
shaded. 
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Table 1 . Globally averaged, area weighted annual mean age of air (years) derived from SF 6 (relative to the 
global mean SF 6 at the ground), at 19, 27, and 35 km, along with the stratospheric steady state lifetimes 
(years) of CFCI 3 and N 2 0, for each transport scenario as discussed in the text. Also shown are the lifetimes, 
along with the 1 standard deviation uncertainties, estimated from ER-2 aircraft data by Volk et al. [1997]. We 
have also included a rough approximation of a global mean age of air, along with a range, derived from OMS 
balloon SF& data. These values are area weighted averages of data taken at 7°S, 35°N, and 65°N. The averages 
and ranges correspond to observations within ±0.5 km of the altitude indicated. 


Mean Age (yrs) Lifetime (yrs) 


Scenario 

Description 

19 km 

27 km 

35 km 

CFCI 3 

N 2 0 

A 

Base (1999 model) 

3.1 

4.4 

4.9 

68 

136 

B 

very weak trop/strat circulation 

5.3 

7.7 

8.4 

108 

195 

C 

weak trop/strat circulation 

4.3 

6.1 

6.7 

89 

166 

D 

strong trop/strat circulation 

2.3 

3.3 

3.9 

55 

118 

E 

very strong trop/strat circulation 

1.5 

2.2 

2.7 

40 

98 

F 

stratospheric minimum K zz — 001 m 2 /s 

3.5 

4.8 

5.3 

71 

139 

G 

stratospheric minimum K zz = .1 m 2 /s 

2.6 

3.9 

4.4 

62 

130 

H 

stratospheric minimum K zz = 1 m 2 /s 

1.0 

1.8 

2.3 

34 

96 

I 

weak mesospheric gravity wave drag, K zz 

3.2 

4.5 

5.1 

69 

142 

J 

strong mesospheric gravity wave drag, K zz 

2.9 

4.1 

4.5 

63 

121 

K 
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